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Abstract 


A search for single top-quark production via flavour-changing neutral current processes from 
gluon plus up- or charm-quark initial states in proton-proton collisions at the LHC is presen¬ 
ted. Data collected with the ATLAS detector in 2012 at a centre-of-mass energy of 8 TeV 
and corresponding to an integrated luminosity of 20.3 fb“^ are used. Candidate events for a 
top quark decaying into a lepton, a neutrino and a jet are selected and classified into signal- 
and background-like candidates using a neural network. No signal is observed and an upper 
limit on the production cross-section multiplied by the t Wb branching fraction is set. 
The observed 95 % CL limit is CTqg^t X 3{t Wb) < 3.4 pb and the expected 95 % CL 
limit is (Tqg^t X !B{t Wb) < 2.9 pb. The observed limit can be interpreted as upper limits 
on tbe coupling constants of the flavour-changing neutral current interactions divided by the 
scale of new physics Kugt!^ < 5.8 X 10 ^TeV * and KcqtIA < 13 X 10 ^TeV * and on the 
branching fractions !B{t ug) < 4.0 x 10”^ and !B{t eg) < 20 x 10”^. 
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1 Introduction 


The top quark is the most massive elementary particle known, with a mass mtop = 173.3 ± 0.8 GeV [1] 
close to the electroweak symmetry breaking scale. This makes it an excellent object with which to test 
the Standard Model (SM) of particle physics, as well as to search for phenomena beyond the SM. 

At the LHC, top quarks are primarily produced in pairs via the strong interaction. In addition to the 
predominant pair-production process, top quarks are produced singly through three different subprocesses 
via the weak interaction: the t-channel, which is the dominant process, involving the exchange of a space¬ 
like W boson; the Wt associated production, involving the production of a real W boson; and the ^-channel 
process involving the production of a time-like W boson. 

As a consequence of the large value, which is close to one, of the Vtb element in the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix, the predominant decay channel of top quarks is t ^ Wb. Transitions between 
top quarks and other quark flavours mediated by neutral gauge bosons, so-called flavour-changing neutral 
currents (FCNC), are forbidden at tree level and suppressed at higher orders in the SM [2]. However, 
several extensions to the SM exist that significantly enhance the production rate and hence the branching 
fractions, S, of FCNC processes. Examples of such extensions are the quark-singlet model [3-5], two- 
Higgs-doublet models with or without flavour conservation [6-11], the minimal supersymmetric standard 
model [12-18] or supersymmetry with /^-parity violation [19, 20], models with extra quarks [21-23], or 
the topcolour-assisted technicolour model [24]. Reviews can be found in Refs. [25] and [26]. Many of 
these models allow for enhanced FCNC production rates, e.g. by permitting FCNC interactions at tree 
level or introducing new particles in higher-order loop diagrams. The predicted branching fractions for 
top quarks decaying to a quark and a neutral boson can be as large as 10“^ to 10“^ for certain regions of 
the parameter space in the models mentioned. However, the experimental limits have not excluded any 
specific extension of the SM for the process t ^ qg so far. 

Among FCNC top-quark decays of the form t ^ qX with X = Z,H, y, g, modes involving a Z boson, a 
Higgs boson (H), or a photon (y) are usually studied directly by searching for final states containing the 
corresponding decay particles. However, the mode t —> qg, where q denotes either an up quark, u, or a 
charm quark, c, is nearly indistinguishable from the overwhelming background of multi-jet production 
via quantum chromodynamic (QCD) processes. For the t ^ qg mode, much better sensitivity can be 
achieved by searching for anomalous single top-quark production {qg —> t) where a u- or c-quark and a 
gluon g, originating from the colliding protons, interact to produce a single top quark. A leading-order 
diagram for top-quark production in the qg ^ t mode as well as a SM decay of the top quark is shown in 
Fig. 1.1 

Anomalous FCNC couplings can be described in a model-independent manner using an effective operator 
formalism [27], which assumes the SM to be the low-energy limit of a more general theory that is valid at 
very high energies. The effects of this theory below a lower energy scale. A, are perceived through a set of 
effective operators of dimension higher than four. The formalism therefore allows the new physics to be 
described by an effective Lagrangian consisting of the SM Lagrangian and a series of higher-dimension 
operators, which are suppressed by powers of 1/A. The new physics scale. A, has a dimension of energy 
and is related to the mass cut-off scale above which fhe effective fheory breaks down, hence characferising 
fhe energy scale at which the new physics manifests itself in the theory. A further method for simplifying 
the formalism is to only consider operators of interest that have no sizeable impact on physics below the 
TeV scale, following Ref. [28]. 

* Charge conjugate production and decay modes are implied throughout this paper. 
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Figure 1; Leading-order Feynman diagram for FCNC top-quark production in the qg ^ t mode followed by the 
decay of the top quark into a h-quark and a W boson, where the W boson decays into a lepton and a neutrino. 


The interest of this paper lies in effective dimension-six operators, which contribute to flavour-changing 
interactions in the strong sector; thus no operators with electroweak gauge bosons are considered. In par¬ 
ticular, the operators describing FCNC couplings to a single top quark are of interest here; they describe 
strong FCNC vertices of the form qgt and can be written as [29] : 

where stands for a right-handed quark singlet, for a left-handed quark doublet, is the gluon 
field strength tensor, <1) the charge conjugate of the Higgs doublet, T" are the Gell-Mann matrices and 
cr^’' is the anti-symmetric tensor. The indices (/, j) of the spinors are flavour indices indicating the quark 
generation. By requiring a single top quark in the interaction, one of the indices can always be set equal 
to 3 while the other index is either 1 or 2. Hence, the remaining fermion field in the interaction is either a 
u- or a c-quark. Apart from direct single top-quark production, these operators give rise to interactions of 
the form gg tq and gq —> tg. The processes considered are a subset of these, where a n-quark, c-quark 
or gluon originating from the colliding protons interacts through an s-, t- or n-channel process to produce 
a single top quark, either via a (2 —> 2) process or without the associated production of additional gluons 
or light quarks via a (2 ^ 1) process. 

The corresponding strong FCNC Lagrangian usually is written as [29] : 

q=M,c 

with the real and positive parameters Kgqt (q = u, c) that relate the strength of the new couplings to the 
strong coupling strength, g^, and where t denotes the top-quark field. The parameters fq and hq are 
real, vector and axial chiral parameters, respectively, which satisfy the relation |/^p -i- \hq^ = 1. This 
Lagrangian contributes to both the production and decay of top quarks. 

Experimental limits on the branching fractions of the FCNC top-quark decay channels have been set by 
experiments at the LEP, HERA, Tevatron and EHC accelerators. At present the most stringent upper limits 
at 95 % confidence level (CE) for the coupling constants Kyqt and Kqgt are Kyqtlrritop < 0.12GeV“* [30] 
(ZEUS, HERA) and S{t qg) < 5.7 x lO’^ (ugt) and S{t qg) < 2.7 x 10“^ (cgt) [31] (ATEAS, 
EHC). In the case of t ^ qZ, upper limits on the branching fractions of the top-quark decay have been 
determined to be qZ) < 0.05 % [32] (CMS, EHC). Einally, the most stringent limit for the decay 

t^qH is measured to be !B{t qH) < 0.79 % [33] (ATEAS, EHC). 
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In the allowed region of parameter space for KggtIA, the FCNC production cross-section for single top 
quarks is of the order of picobarns, while the branching fraction for FCNC decays is very small, i.e. 
below 1 %. Top quarks are therefore reconstructed in the SM decay mode t Wb. The W boson can 
decay into a quark-antiquark pair {W —> qiqi) or a charged lepton-neutrino pair (W —> iv)', only the 
latter is considered here. This search targets the signature from the qg ^ t ^ W{-^ €v) b process. 
Events are characterised by an isolated high-energy charged lepton (electron or muon), missing trans¬ 
verse momentum from the neutrino and exactly one jet produced by the hadronisation of the Z^-quark. 
Events with a W boson decaying into a r lepton, where the r decays into an electron or a muon, are also 
included. Several SM processes have the same final-state topology and are considered as background to 
the ECNC analysis. The main backgrounds are F-rJets production (especially in association with heavy 
quarks), where V denotes a IT or a Z boson, SM top-quark production, diboson production, and multi-jet 
production via QCD processes. The studied process can be differentiated from SM single top-quark pro¬ 
duction, which is usually accompanied by additional jets. Eurthermore, ECNC production has kinematic 
differences from the background processes, such as lower transverse momenta of the top quark. 

This paper is organised as follows: Section 2 provides a description of the ATEAS detector. Section 3 
gives an overview of the data and Monte Carlo (MC) samples used for the simulation of signal and 
expected background events from SM processes. In Sect. 4 the event selection is presented. The methods 
of event classification into signal- and background-like events using a neural network are discussed in 
Sect. 5 and sources of systematic uncertainty are summarised in Sect. 6 . The results are presented in 
Sect. 7 and the conclusions are given in Sect. 8 . 


2 ATLAS detector 

The ATEAS detector [34] is a multipurpose collider detector built from a set of sub-detectors, which cover 
almost the full solid angle around the interaction point.^ It is composed of an inner tracking detector (ID) 
close to the interaction point surrounded by a superconducting solenoid providing a 2 T axial magnetic 
field, electromagnetic and hadronic calorimeters, and a muon spectrometer (MS). The ID consists of a 
silicon pixel detector, a silicon microstrip detector providing tracking information within pseudorapidity 
\q\ < 2.5, and a straw-tube transition radiation tracker that covers I 77 I < 2.0. The central electromag¬ 
netic calorimeter is a lead and liquid-argon (EAr) sampling calorimeter with high granularity, and is 
divided into a barrel region that covers \q\ < 1.475 and endcap regions that cover 1.375 < \q\ < 3.2. 
An iron/scintillator tile calorimeter provides hadronic energy measurements in the central pseudorapidity 
range. The endcap and forward regions are instrumented with EAr calorimeters for both the electromag¬ 
netic and hadronic energy measurements, and extend the coverage to \q\ - 4.9. The MS covers \q\ < 2.7 
and consists of three large superconducting toroids with eight coils each, a system of trigger chambers, 
and precision tracking chambers. 


^ ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point in the centre of the detector and 
the z-axis is along the beam direction; the x-axis points towards the centre of the LHC ring and the ^-axis points upwards. 
The pseudorapidity r] is defined as 77 = - ln[tan(6/2)], where the polar angle 6 is measured with respect to the z-axis. The 
azimuthal angle, 0, is measured with respect to the x-axis. Transverse momentum and energy are defined as pj = p sin d and 
Ej = E sin 9, respecfively. The AR distance in ( 77 ,space is defined as AR = ^(Arf)^ 
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3 Data and simulated samples 


This analysis is performed using ^/s = S TeV proton-proton (pp) collision data recorded by the ATLAS 
experiment in 2012. Stringent detector and data quality requirements are applied, resulting in a data 
sample with a total integrated luminosity of 20.3 fb“^ 


3.1 Trigger requirements 

ATLAS employs a three-level trigger system for selecting events to be recorded. The first level (LI) is 
built from custom-made hardware, while the second and third levels are software based and collectively 
referred to as the high-level trigger (HLT). The datasets used in this analysis are defined by high-pT 
single-elecfron or single-muon triggers [35, 36]. 

For the LI calorimeter trigger, which is based on reduced calorimetric information, a cluster in the elec¬ 
tromagnetic calorimeter is required with Ej > 30GeV or with Ej > ISGeV. The energy deposit must 
be well separated from other clusters. At the HLT, the full granularity of the calorimeter and tracking 
information is available. The calorimeter cluster is matched to a track and the trigger electron candidate 
is required to have Ej- > 60 GeV or Ej > 24 GeV with additional isolation requirements. 

The single-muon trigger is based on muon candidates reconstructed in the MS. The triggered events 
require a LI muon trigger-chamber track with a 15 GeV threshold on the pj of the track. At the HLT, the 
requirement is tightened to p^ > 24 GeV with, or 36 GeV without, an isolation criterion. 

3.2 Simulated events 

Simulated event samples are used to evaluate signal and background efficiencies and uncertainties as well 
as to model signal and background shapes. 

For the direct production of top quarks via FCNC, MEtop [29] is used for simulating strong FCNC pro¬ 
cesses at next-to-leading order (NLO) in QCD. It introduces strong top-quark FCNC interactions through 
effective operators. By comparing kinematic distributions for different FCNC couplings, it has been veri¬ 
fied fhaf fhe kinematics of the signal process are independent of the a priori unknown FCNC coupling 
strength. As a conservative approach, only left-handed top quarks (as in the SM) are produced, and the 
decay of the top quark is assumed also to be as in the SM.^ The CTIO [37] parton distribution func¬ 
tion (PDF) sets are used for the generation of the signal events and the renormalisation and factorisation 
scales are set to the top-quark mass. 

The PowHEG-Box [38] generator with the CTIO PDF sets is used to generate it [39] and electroweak single 
top-quark production in the t-channel [40], ^-channel [41] and IVt-channel [42]. All processes involving 
top quarks, including the strong FCNC processes, are produced assuming mtop = 172.5 GeV. The parton 
shower and the underlying event are added using Pythia 6.426 [43], where the parameters controlling the 
modelling are set to the values of the Perugia 2011C tune [44]. 

Vector-boson production in association with jets (V-i-jets) is simulated using the multi-leg leading-order 
(LO) generator Sherpa 1.4.1 [45] with its own parameter tune and the CTIO PDF sets. Sherpa is used 


^ A right-handed top quark would give rise to different angular distributions and hence be easier to separate from SM production. 
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not only to generate the hard process, but also for the parton shower and the modelling of the under¬ 
lying event. IV-f-jets and Z-f-jets events with up to five additional partons are generated. The CKKW 
method [46] is used to remove overlap between partonic configurations generated by the matrix element 
and by parton shower evolution. Double counting between the inclusive V+n parton samples and samples 
with associated heavy-quark pair production is avoided consistently by using massive c- and f7-quarks in 
the shower. 

Diboson events {WW, WZ and ZZ) are produced using Alpgen 2.14 [47] and the CTEQ6L1 PDF sets [48]. 
The partonic events are showered with Herwig 6.5.20 [49], and the underlying event is simulated with 
the Jimmy 4.31 [50] model using the ATLAS Underlying Event Tune 2 [51]. 

All the generated samples are passed through the full simulation of the ATLAS detector [52] based on 
Geant4 [53] and are then reconstructed using the same procedure as for data. The simulation includes 
the effect of multiple pp collisions per bunch crossing. The events are weighted such that the average 
distribution of the number of collisions per bunch crossing is the same as in data. In addition, scale 
factors are applied to the simulated events to take into account small differences observed befween fhe 
efficiencies for fhe frigger, lepfon idenfificafion and fi-quark jef idenfificafion. These scale factors are 
defermined using confrol samples. 


4 Event selection 

The expecfed signafure of signal evenfs is used fo perform fhe even! selecfion. Evenfs confaining ex- 
acfly one isolated elecfron or muon, missing fransverse momenfum and one jef, which is required to be 
idenlified as a jef originafing from a fi-quark, are selected. 


4.1 Object definition and event selection 

Elecfron candidates are selected from energy deposifs (clusfers) in fhe LAr elecfromagnefic calorimefer 
associafed wifh a well-measured frack fulfilling sfricf qualify requiremenfs [54]. Elecfron candidates are 
required fo satisfy px > 25 GeV and |?7ciusl < 2.47, where pdus denofes fhe pseudorapidify of fhe cluster. 
Clusfers falling in fhe calorimefer barrel-endcap fransifion region, corresponding fo 1.37 < liydusl < 1-52, 
are ignored. High-px elecfrons associafed wifh fhe VT-boson decay can be mimicked by hadronic jefs 
reconsfrucfed as elecfrons, elecfrons from fhe decay of heavy quarks, and phofon conversions. Since 
elecfrons from fhe IT-boson decay are fypically isolafed from hadronic jef acfivify, backgrounds can be 
suppressed by isolation criferia, which require minimal calorimefer acfivify and only allow low-px fracks 
in an cone around fhe election candidafe. Isolafion cufs are opfimised fo achieve a uniform cuf 
efficiency of 90 % as a function of pcius and fransverse energy, Ex- The direction of fhe elecfron candidafe 
is faken as fhaf of fhe associafed frack. For fhe calorimefer isolafion a cone size of AT? = 0.2 is used. 
In addition, fhe scalar sum of all frack fransverse momenfa wifhin a cone of size A/? = 0.3 around fhe 
elecfron direcfion is required fo be below a px-dependenf fhreshold in fhe range befween 0.9 GeV and 
2.5 GeV. The frack belonging fo fhe election candidafe is excluded from Ibis requiremenf. 

Muon candidafes are reconsfrucfed by mafching frack segmenfs or complete fracks in fhe MS wifh fracks 
found in fhe ID [55]. The final candidafes are required fo have a fransverse momenfum px > 25 GeV and 
fo be in fhe pseudorapidify region |p| < 2.5. Isolafion criferia are applied to reduce background evenfs in 
which a high-px muon is produced in fhe decay of a heavy-flavour quark. An isolafion variable [56] is 
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defined as the scalar sum of the transverse momenta of all tracks with pj above 1 GeV, except the one 
matched to the muon, within a cone of size A/?iso - lOG&V/pjijS). Muon candidates are accepted if they 
have an isolation to pT:{p) ratio of less than 0.05. An overlap removal is applied between the electrons 
and the muons, rejecting the event if the electron and the muon share the same ID track. 

Jets are reconstructed using the anti-k, algorithm [57] with a radius parameter of 0.4, using topolo¬ 
gical clusters [58] as inputs to the jet finding. The clusters are calibrated with a local cluster weighting 
method [59]. Calibrated jets using an energy- and //-dependent simulation-based calibration scheme, with 
in situ corrections based on data, are at first required to have pj > 25 GeV and \r]\ < 2.5. The jet energy 
is further corrected for the effect of multiple pp interactions, both in data and in simulated events. 

If any jet is within AR = 0.2 of an electron, the closest jet is removed, since in these cases the jet and 
the electron are very likely to correspond to the same physics object. Remaining electron candidates 
overlapping with jets within a distance AR < 0.4 are subsequently rejected. To reject jets from pile-up 
events, a so-called jet-vertex fraction criterion is applied for jets with pj < 50 GeV and |? 7 | < 2.4: at least 
50 % of the scalar sum of the pj of tracks within a jet is required to be from tracks compatible with the 
primary vertex^ associated with the hard-scattering collision. The final selected jet is required to have 
Pi > 30 GeV and must also be identified as a jet originating from a t/-quark (^/-tagged). 

In this analysis, a ^/-tagging algorithm that is optimised to improve the rejection of c-quark jets is used, 
since W+c production is a major background. A neural-network-based algorithm is used, which combines 
three different algorithms exploiting the properties of a t/-hadron decay in a jet [60]. The chosen working 
point corresponds to a ^/-tagging ethciency of 50 %, when cutting on the discriminant, and a c-quark jet 
and light-parton jet mistag acceptance of 3.9 % and 0.07 %, respectively, as measured in tt events [61, 
62]. 

The missing transverse momentum (with magnitude S™®®) is calculated based on the vector sum of energy 
deposits in the calorimeter projected onto the transverse plane [63]. All cluster energies are corrected 
using the local cluster calibration scheme. Clusters associated with a high-pj jet or electron are further 
calibrated using their respective energy corrections. In addition, contributions from the pi of selected 
muons are included in the calculation of Due to the presence of a neutrino in the final state of the 
signal process, > 30 GeV is required. Lepton candidates in multi-jet events typically arise from 
charged tracks being misidentified as leptons, electrons arising from converted photons and leptons from 
c- and (/-hadron decays. Such candidates are collectively referred to as fake leptons. As such, the multi¬ 
jet events tend to have low L"™®® and low IT-boson transverse mass,^ mi{W), relative to single top-quark 
events. Therefore, an additional requirement on mi{W) is an effective way to reduce this background. 
The selection applied is mi{W) > 50 GeV. In order to further suppress the multi-jet background and also 
to remove poorly reconstructed leptons with low transverse momentum, a requirement on the transverse 
momentum of leptons and the azimuthal angle between the lepton and jet is applied: 


p\ > 90 GeV I 


n-- |A0(Ajet)| 
n -1 


( 1 ) 


The parameters of the cut are motivated by the distribution of multi-jet events, obtained in the signal 
region, where the simulated backgrounds except the multi-jet contribution are subtracted from data. Al¬ 
most no signal events are removed by this cut. The distribution of the transverse momentum of the lepton 
versus the azimuthal angle between the lepton and the jet is shown in Fig. 2. 

^ The primary vertex is defined as the vertex with the largest 2 //| of the associated tracks. 

^ The VT-boson transverse mass is defined as: m^iW) = .^2 - Pt(£) ■ where Pt(C) denotes the transverse 
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Figure 2: The transverse momentum of the lepton versus the azimuthal angle between the lepton and the jet. The 
colours indicate the number of events in data after the simulated backgrounds except the multi-jet contribution have 
been subtracted and before the cut given by Eq. 1 is applied. The solid black line shows the cut. 


In addition to the signal region defined by this selection, a control region is defined wifh fhe same kin¬ 
ematic requiremenfs, buf wifh a less sfringenf fj-fagging requiremenf wifh an efficiency of 85 %, and 
excluding evenfs passing fhe fighfer signal-region fj-fagging selection. This confrol region is designed 
such fhaf fhe resulfing sample is dominafed by IT-f-jels producfion, which is fhe dominanf background. 


4.2 Background estimation 

For all background processes excepf fhe mulfi-jef background, fhe normalisafions are esfimafed by using 
Monfe Carlo simulafion scaled fo fhe fheorefical cross-secfion predicfions, using rriiop = 172.5 GeV. In 
order fo check fhe modelling of kinemafic disfribufions, correction factors fo fhe normalisation of fhe 
VF-i-jefs and tt and single-fop processes are subsequenfly defermined simulfaneously in fhe confexf of fhe 
mulfi-jef background estimation. 

The SM single fop-quark producfion cross-sections are calculafed fo approximafe nexf-fo-nexf-fo-leading- 
order (NNLO) precision. The producfion via fhe f-channel exchange of a virfual W boson has a predicfed 
cross-secfion of 87 pb [64]. The cross-secfion for fhe associafed producfion of an on-shell W boson and 
a fop quark (Wt channel) has a predicted value of 22.3pb [65], while fhe 5-channel producfion has a 
predicfed cross-secfion of 5.6pb [66]. The resulfing weighfed average of fhe fheorefical uncerfainfies 
including PDF and scale uncerfainfies of fhese fhree processes is 10 %. 

The cross-section of fhe tt process is normalised fo 238 pb, calculafed af NNLO in QCD including re- 
summafion of nexf-lo-next-lo-leading logarifhmic (NNLL) soff gluon terms [67-71] wifh Top-t-t2.0 [72]. 
The PDF and as uncerfainfies are calculafed using fhe PDF4LHC prescripfion [73] wifh fhe MSTW2008 


momentum of the lepton and piff) = IfhCDI- 
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NNLO [74, 75] at 68 % CL, the CTIO NNLO [37, 76], and the NNPDF 2.3 [77] PDF sets, and are added 
in quadrature to the scale uncertainty, yielding a final uncertainty of 6 %. 

The cross-sections for inclusive W- and Z-boson production are predicted with NNLO precision using 
the FEWZ program [78, 79], resulting in a LO-to-NNLO .K'-factor of 1.10 and an uncertainty of 4 %. The 
uncertainty includes the uncertainty on the PDF and scale variations. The scale factor is applied to the 
prediction based on the LO Sherpa calculation and the flavour composifion is also faken from fhe MC 
samples. The modelling of fhe fransverse momenfum of fhe W boson in fhe fT-i-jefs sample is improved 
by reweighfing fhe simulafed samples fo dafa in fhe IT-i-jefs-dominafed confrol region. 

LO-fo-NLO A'-facfors obfained wifh MCFM [80] of fhe order of 1.3 are applied fo fhe Alpgen LO pre¬ 
dictions for diboson producfion. Since fhe diboson process is freafed fogefher wifh Z-boson producfion in 
fhe sfafisfical analysis and fhe fraction of selecfed evenfs is only 5 %, fhe same uncerfainfies as used for 
fhe Z-i-jefs process are assumed. 

Mulfi-jef evenfs may be selecfed if a jef is misidenfified as an isolafed lepfon or if fhe evenf has a non- 
prompf lepfon fhaf appears fo be isolafed. The normalisafion of fhis background is obfained from a fif fo 
fhe observed F™'®® disfribufion, performed bofh in fhe signal and confrol regions. In order fo consfrucf a 
sample of mulfi-jef background evenfs, differenf mefhods are adopfed for fhe elecfron and muon channels. 
The ‘jef-lepfon’ model is used in fhe elecfron channel while fhe ‘anfi-muon’ model is used in fhe muon 
channel [81]. In fhe jef-lepfon model, a shape for fhe mulfi-jef background is esfablished using evenfs 
from a Pythia dijef sample, which are selecfed using same criferia as fhe sfandard selecfion, buf wifh a jef 
used in place of fhe elecfron candidafe. Each candidafe jef has fo fulfil fhe same pj and rj requiremenfs 
as a sfandard lepfon and deposif 80-95 % of ifs energy in fhe elecfromagnefic calorimefer. Evenfs wifh an 
elecfron candidafe passing fhe elecfron cufs described in Secf. 4.1 are rejecfed and an evenf is accepfed if 
exacfly one ‘jef-lepfon’ is found. The anfi-muon model is derived from collision dafa. In order fo selecf 
a sample fhaf is highly enriched wifh muons from mulfi-jef evenfs, some of fhe muon identification cufs 
are inverfed or changed, e.g. fhe isolation criferia are inverfed. 

To defermine fhe normalisafion of fhe mulfi-jef background femplafe, a binned maximum-likelihood fif is 
performed on fhe disfribufion using fhe observed dafa, affer applying all selecfion criferia excepf for 
fhe cuf on E™*®. Eifs are performed separafely in fwo p regions for elections: in fhe endcap (|? 7 | > 1.52) 
and cential (|7/| < 1.37) region of fhe elecfromagnefic calorimefer, i.e. fhe tiansifion region is excluded. 
Eor muons, fhe complefe p region is used. The mulfi-jef femplafes for bofh fhe elections and fhe muons 
are fitted fogefher wifh femplafes derived from MC simulation for all ofher background processes (fop 
quark, VT-i-lighf flavour (EE), IT-i-heavy flavour (HE), Z-i-jels, dibosons). Accepfance uncerfainfies are 
accounfed for in fhe tiffing process in fhe form of addifional consfrained nuisance paramefers. Eor fhe 
purpose of fhese fifs, fhe confribufions from IT-i-EE and VT-i-HE, fhe confribufions from ti and single fop- 
quark producfion, and fhe confribufions from Z-i-jefs and diboson producfion are each combined info one 
femplafe. The normalisafion of fhe femplafe for Z-i-jefs and diboson producfion is fixed during fhe fif, as 
ifs confribufion is very small. 

The E™®® disttibufions affer rescaling fhe differenf backgrounds and fhe mulfi-jefs femplafe fo fheir re¬ 
spective fif resulfs are shown in Eig. 3 for bofh fhe elecfron and fhe muon channels. The tiffed scale faclors 
for fhe ofher femplafes are close fo 1. 

Table 1 provides fhe evenf yields affer fhe complefe evenf selecfion for fhe confrol and signal regions. The 
yields are calculafed using fhe accepfance from MC samples normalised fo fheir respective fheorefical 
cross-sections including fhe (N)NEO E-facfors, while fhe number of expecfed evenfs for fhe mulfi-jef 
background is obfained from fhe maximum-likelihood fif. Each evenf yield uncerfainfy combines fhe 
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Figure 3: Fitted distributions of the missing transverse momentum for (a) central electrons and (b) muons in 
the control region and for (c) central electrons and (d) muons in the signal region. The last histogram bin includes 
overflow events and the hatched error bands contain the MC statistical uncertainty combined with the normalisation 
uncertainty on the multi-jet background. 
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statistical uncertainty, originating from the limited size of the simulation samples, with the uncertainty 
on the cross-section or normalisation. The observed event yield in data agrees well with the background 
prediction. For comparison, a 1 pb FCNC cross-section would lead to 530 events in the signal region. 
The corresponding efficiency for selecting FCNC events is 3.1 %. 


Process 

Control region 

Signal region 

Single top 

11500± 620 

14 400 ± 

770 

tt 

10 700 ± 650 

12 000 ± 

740 

W+LF 

526 900 ±130000 

6 700± 

1900 

W+UF 

445 200 ±240 000 

62 100 ±34 000 

Z-i-jets 

40 000 ± 9 700 

4 990± 

1200 

Multi-jet 

68 300 ± 12 000 

7 430± 

1300 

Total expected 

1 100 000 ±280 000 

107 000 ±34 000 

Data 

1 112225 

108152 



Table 1; Number of observed and expected events in the control and signal region for all lepton categories added 
together. The uncertainties shown are derived using the statistical uncertainty from the limited size of the samples 
and the uncertainty on the theoretical cross-section only or multi-jet normalisation. The scale factors obtained from 
the multi-jet background fit are not applied when determining the expected number of events. 

Kinematic distributions in the control region of the identified lepton, reconstructed jet, and mT(lT) 
are shown in Fig. 4 for the combined electron and muon channels. These distributions are normalised 
using the scale factors obtained in the fit to estimate the multi-jet background. Overall, good agree¬ 
ment between the observed and expected distributions is seen. The trends that can be seen in some of the 
distributions are covered by the systematic uncertainties. 


5 Analysis strategy 


As no single variable provides sufficient discrimination between signal and background events and the 
separation power is distributed over many correlated variables, multivariate analysis techniques are ne¬ 
cessary to separate signal candidates from background candidates. A neural-network (NN) classifier [82] 
that combines a three-layer feed-forward neural network with a preprocessing of the input variables is 
used. The network infrastructure consists of one input node for each input variable plus one bias node, an 
arbitrary number of hidden nodes, and one output node, which gives a continuous output in the interval 
[-1,1]. The training is performed with a mixture of 50 % signal and 50 % background events, where the 
different background processes are weighted according to their number of expected events. Only pro¬ 
cesses from simulated events are considered in the training, i.e. no multi-jet events are used. In order to 
check that the neural network is not overtrained, 20 % of the available simulated events are used as a test 
sample. Subsequently, the NN classifier is applied to all samples. 

The —> f —> Mv process is characterised by three main differences from SM processes. Firstly, the pT 
distribution of the top quark is much softer than the pj distribution of top quarks produced through SM 
top-quark production, since the top quark is produced almost without transverse momentum. Hence, the 
W boson and Zj-quark from the top-quark decay are produced almost back-to-back in the transverse plane. 
Secondly, unlike in the IF/Z-i-jets and diboson backgrounds, the W boson from the top-quark decay has a 
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Figure 4: Kinematic distributions in the control region for the combined electron and muon channels. All processes 
are normalised to the result of the binned maximum-likelihood ht used to determine the fraction of multi-jet events. 
Shown are: (a) the transverse momentum and (b) pseudorapidity of the lepton, (c) the transverse momentum and (d) 
pseudorapidity of the jet, (e) the missing transverse momentum and (f) VF-boson transverse mass. The last histogram 
bin includes overflow events and the hatched band indicates the combined statistical and systematic uncertainties, 
evaluated after the ht discussed in Section 7. 
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high momentum and its decay products tend to have small angles. Lastly, the top-quark charge asymmetry 
differs between FCNC processes and SM processes in the ugt channel. In pp collisions, the FCNC 
processes are predicted to produce four times more single top quarks than anti-top quarks, whereas in SM 
single top-quark production and in all other SM backgrounds this ratio is at most two. Several categories 
of variables are considered as potential discriminators between the signal and background processes. 
Apart from basic event kinematics such as the mjCW) or //t (the scalar sum of the transverse momenta 
of all objects in the final sfafe), various objecf combinafions are considered as well. These include fhe 
basic kinematic properfies of reconsfrucfed objecfs like fhe W boson and fhe fop quark, as well as angular 
disfances in 77 and (p befween fhe reconsfrucfed and final-sfafe objecfs in fhe laboratory frame and in fhe 
resf frames of fhe W boson and fhe lop quark. In order fo reconsfrucf fhe four-veclor of fhe W boson, 
a mass consfrainf is used. A defailed descripfion of the top-quark reconstruction is given in Ref. [83]. 
Further, integer variables such as the charge of the lepton are considered. 

The ranking of the variables in terms of their discrimination power is automatically determined as part of 
the preprocessing step and is independent of the training procedure [84].^ Only the highest-ranking vari¬ 
ables are chosen for the training of the neural network. Each variable is tested beforehand for agreement 
between the background model and the distribution of the observed events in the control region. Using 
only variables with an a priori defined separafion power, 13 variables remain in the network. Table 2 
shows a summary of the variables used, ordered by their importance. The probability density of the three 
most important discriminating variables for the dominant background processes together with the signal 
is displayed in Fig. 5. 

The distributions for three of the four most important variables in the control and signal regions are shown 
in Fig. 6 . The shape of the multi-jet background is obtained using the samples described in Sec. 4.2. The 
distribution of is shown in Fig. 4(a) for the control region. The distributions are normalised using the 
scale factors obtained in the binned maximum-likelihood fit to the distribution. 


Variable Definition 


mT(top) 

Ft 

A/?(top, C) 


^-iet 

Ft 

A0(top, b-]Qt) 
cos 6{(, b-]et) 

mi{W) 

77 ^ 

A 0 (top, W) 
A/?(top, f 7 -jet) 

77top 

w 

Ft 


Transverse mass of the reconstructed top quark 
Transverse momentum of the charged lepton 

Distance in the r]-<p plane between the reconstructed top quark and the charged 
lepton 

Transverse momentum of the f 7 -tagged jet 

Difference in azimufh befween fhe reconsfrucfed fop quark and fhe f 7 -tagged jef 
Opening angle of the three-vectors between the charged lepton and the f 7 -tagged jet 
Charge of the lepton 
IT-boson transverse mass 
Pseudorapidity of the charged lepton 

Difference in azimufh befween fhe reconsfrucfed fop quark and fhe W boson 
Disfance in the 77 -^ plane between the reconstructed top quark and the f 7 -tagged jet 
Pseudorapidity of the reconstructed top quark 
Transverse momentum of the W boson 


Table 2: Variables used in the training of the neural network ordered by their descending importance. 


® The ranking is done according to the correlation to the output. 
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Figure 5: Probability densities of the three most important discriminating variables; (a) the transverse mass of the 
reconstructed top quark; (b) the transverse momentum of the charged lepton; and (c) the distance in the ri-([> plane 
between the charged lepton and the reconstructed top quark. The last histogram bin includes overflows. 


The resulting neural-network output distributions for the most important background processes and the 
signal are displayed in Fig. 7 as probability densities and in Figs. 8(a) and 8(b) normalised to the number 
of expected events in the control and signal regions, respectively. Signal-like events have output values 
close to 1, whereas background-like events accumulate near -1. Overall, good agreement within system¬ 
atic uncertainties between data and the background processes is observed in both the control and signal 
regions. 


6 Systematic uncertainties 


Systematic uncertainties are assigned to account for detector calibration and resolution uncertainties, as 
well as the uncertainties on theoretical predictions. These can affect the normalisation of the individual 
backgrounds and the signal acceptance (acceptance uncertainties) as well as the shape of the neural- 
network output distribution (shape uncertainties). Quoted relative uncertainties refer to acceptance of the 
respective processes unless stated otherwise. 


6.1 Object modelling 

The effects of the systematic uncertainties due to the residual differences between data and Monte Carlo 
simulation, uncertainties on jets, electron and muon reconstruction after calibration, and uncertainties on 
scale factors that are applied to the simulation are estimated using pseudo-experiments. 

Uncertainties on the muon (electron) trigger, reconstruction and selection efficiency scale factors are 
estimated in measurements of Z ^ ///r (Z —> ee and W —> ev) production. The scale factor uncertainties 
are as large as 5 %. To evaluate uncertainties on the lepton momentum scale and resolution, the same 
processes are used [85]. The uncertainty on the charge misidentification acceptances were studied and 
found to be negligible for this analysis. 

The jet energy scale (JES) is derived using information from test-beam data, LHC collision data and 
simulation. Its uncertainty varies between 2.5 % and 8 %, depending on jet pj and p [59]. This includes 
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Figure 6; Distributions of three important discriminating variables (except for the transverse momentum of the 
lepton); (a) and (d) the top-quark transverse mass in the control and signal regions; (b) and (e) the AR between 
the lepton and the reconstructed top quark in the control and signal regions; (c) and (f) the Acp between the jet and 
the reconstructed top quark. All processes are normalised using the scale factors obtained in the binned maximum- 
likelihood ht to the distribution. The FCNC signal cross-section is scaled to 50 pb and overlayed on the 
distributions in the signal region. The last histogram bin includes overflow events and the hatched band indicates 
the combined statistical and systematic uncertainties, evaluated after the fit discussed in Section 7. 
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Figure 7; Probability density of the neural-network output distribution for the signal and the most important back¬ 
ground processes. 
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Figure 8: Neural-network output distribution (a) in the control region and (b) in the signal region. The shape of the 
signal scaled to 50 pb is shown in (b). All background processes are shown normalised to the result of the binned 
maximum-likelihood ht used to determine the fraction of multi-jet events. The hatched band indicates the combined 
statistical and systematic uncertainties, evaluated after the ht discussed in Section 7. 
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uncertainties in the fraction of jets induced by gluons and mismeasurements due to close-by jets. Addi¬ 
tional uncertainties due to pile-up can be as large as 5 %. An additional jet energy scale uncertainty of up 
to 2.5 %, depending on the pj of the jet, is applied for ft-quark-induced jets due to differences between 
light-quark and gluon jets compared to jets containing Z^-hadrons. Additional uncertainties are from the 
modelling of the jet energy resolution and the missing transverse momentum, which accounts for con¬ 
tributions of calorimeter cells not matched to any jets, soft jets, and pile-up. The effect of uncertainties 
associated with the jet-vertex fraction is also considered for each jet. 

Since the analysis makes use of Z^-tagging, the uncertainties on the b- and c-tagging efficiencies and the 
mistag acceptance [61, 62] are taken into account. 


6.2 Multi-jet background 


For the multi-jet background, an uncertainty on the estimated multi-jet fractions and the modelling is 
included. The systematic uncertainty on the fractions, as well as a shape uncertainty, are obtained by 
comparing to an alternative method, the matrix method [81]. The method estimates the number of multi¬ 
jet background events in the signal region based on loose and tight lepton isolation definitions, the latter 
selection being a subset of the former. The number of multi-jet events passing the tight (signal) 
isolation requirements can be expressed as: 


^flight _ Cfake 

Creal “ ffake 




where treai and efake are the efficiencies for real and fake loose leptons being selected as tight leptons, 
^loose number of selected events in the loose sample, and is the number of selected events 
in the signal sample. By comparing the two methods, the uncertainty on the fraction of multi-jet events 
is estimated to be 17 %. The shape uncertainty is constructed by comparing the neural-network output 
distributions of the jet-lepton and anti-muon samples with the distributions obtained using the matrix 
method. 


6.3 Monte Carlo generators 

Systematic effects from the modelling of the signal and background processes are taken into account by 
comparing different generator models and varying the parameters of the event generation. The effect 
of parton-shower modelling for the top-quark processes is tested by comparing two Powheg samples 
interfaced to Herwig and Pythia, respectively. There are also differences associated with the way in 
which double-counted events in the NLO corrections and the parton showers are removed. These are 
estimated by comparing samples produced with the MC@NLO method and the Powheg method. 

The difference between the top-quark mass used in the simulations and the measured value has negligible 
effect on the results. 

For the single top-quark processes, variations of initial- and final-state radiation (ISR and FSR) together 
with variations of the hard-process scale are studied. The uncertainty is estimated using events generated 
with Powheg interfaced to Pythia. Factorisation and renormalisation scales are varied independently by 
factors of 0.5 and 2.0, while the scale of the parton shower is varied consistently with the renormalisation 
scale using specialised Perugia 2012 tunes [44]. The uncertainty on the amounts of ISR and FSR in the 
simulated tt sample is assessed using Alpgen samples, showered with Pythia, with varied amounts of 
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initial- and final-state radiation, which are compatible with the measurements of additional jet activity in 
tt events [86]. 

The effect of applying the IT-boson px reweighting was studied and found to have negligible impact on 
the shape of the neural-network output distribution and the measured cross-section. Hence no systematic 
uncertainty due to this was assigned. 

Finally, an uncertainty is included to account for statistical effects from the limited size of the MC 
samples. 


6.4 Parton distribution functions 

Systematic uncertainties related to the parton distribution functions are taken into account for all samples 
using simulated events. The events are reweighted according to each of the PDF uncertainty eigen¬ 
vectors or replicas and the uncertainty is calculated following the recommendation of the respective PDF 
group [73]. The final PDF uncerfainfy is given by fhe envelope of fhe esfimafed uncerfainfies for fhe CTIO 
PDF sef, fhe MSTW2008 PDF sef and fhe NNPDF 2.3 PDF sef. 

6.5 Theoretical cross-section normalisation 

The fheorefical cross-secfions and fheir uncerfainfies are given in Secf. 4.2 for each background process. 
Since fhe single fop-quark t-, Wt-, and 5-channel processes are grouped fogefher in fhe sfafisfical analysis, 
fheir uncerfainfies are added in proportion fo fheir relafive fracfions, leading fo a combined uncerlainfy of 
10 %. 

A cross-secfion uncerfainfy of 4 % is assigned for fhe W/Z-i-fO jef) process, while ALPGEN paramefer 
variations of fhe faclorisalion and renormalisafion scale and fhe mafching paramefer consisfent wifh ex¬ 
perimental data yield an uncertainty on the cross-section ratio of 24 %. For IT-i-HF production, a conser¬ 
vatively estimated uncertainty on the HF fraction of 50 % is added. This uncertainty is also applied to the 
combined Z-i-jets and diboson background. 


6.6 Luminosity 

The uncertainty on the measured luminosity is estimated to be 2.8 %. It is derived from beam-separation 
scans performed in November 2012, following the same methodology as that detailed in Ref. [87]. 


7 Results 

In order to estimate the signal content of the selected sample, a binned maximum-likelihood fit to the 
complete neural-network output distributions in the signal region is performed. Including all bins of the 
neural-network output distributions in the fit has the advantage of making maximal use of all signal events 
remaining after the event selection, and, in addition, allows the background acceptances to be constrained 
by the data. 
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The signal rates, the rate of the single top-quark and tt background and the rate of the IT-i-HF background 
are fitted simultaneously. The event yields of the multi-jet background, the IT-i-LF and the combined 
Z-i-jets/diboson background are not allowed to vary in the fit, but instead are fixed fo fhe esfimafes given 
in Table 1 . 

No significanl rale of FCNC single lop-quark producfion is observed. An upper limil is sel using hy- 
polhesis lesfs. The compalibilily of fhe dala wilh fhe signal hypolhesis, which depends on fhe coupling 
consfanfs, and fhe background hypolhesis is evalualed by performing a frequenlisl hypolhesis lesl based 
on pseudo-experimenls, corresponding lo an integrated luminosity of 20.3 fb“^ Two hypolheses are 
compared: Ihe null hypolhesis, Hq, and Ihe signal hypolhesis, H\, which includes FCNC single lop-quark 
producfion. For bolh scenarios, ensemble lesls, i.e. large sels of pseudo-experimenls, are performed. Sys- 
lemalic uncerlainlies are included in Ihe pseudo-experimenls using variations of Ihe signal acceplance, 
Ihe background acceplances and Ihe shape of Ihe neural-nelwork oulpul dislribulion due lo all sources of 
uncertainly described in Ihe previous section. 

To distinguish belween Ihe Iwo hypolheses, Ihe so-called Q value is used as a lesl slalislic. If is defined 
as Ihe ratio of Ihe likelihood function L, evaluated for Ihe differenl hypolheses: 



( 2 ) 


)FCNC 


where yS 


is Ihe scale factor for Ihe number of evenls expected from Ihe signal process for an assumed 


production cross-section. Systematic uncerlainlies are included by varying Ihe predicted number of evenls 
for Ihe signal and all background processes in Ihe pseudo-experimenls. 

The CL^ melhod [88] is used to derive confidence levels (CL) for a cerlain value of and A 

particular signal hypolhesis Hi, determined by given coupling conslanls KugtIA and Kcgt!A, is excluded 
al Ihe 95 % CL if a CL^ < 0.05 is found. The observed 95 % CL upper limil on Ihe anomalous FCNC 
single lop-quark production cross-section multiplied by Ihe t —> Wb branching fraction, including all 
uncerlainlies, is 3.4 pb, while Ihe expected upper limil is 2.9]|;j-2 pb. 

To visualise Ihe observed upper limil in Ihe neural-nelwork oulpul dislribulion, Ihe FCNC signal process 
scaled to 3.4 pb slacked on top of all background processes is shown in Fig. 9. 

The lolal uncerlainly is dominated by Ihe jel energy resolution uncertainly, Ihe modelling of and 
Ihe uncertainty on Ihe normalisation and Ihe modelling of Ihe mulli-jel background. A summary of all 
considered sources and Iheir impacl on Ihe expected upper limil is shown in Table 3. 

Using Ihe NLO predictions for Ihe FCNC single lop-quark production cross-section [89,90] and assuming 
!B{t —> Wb) - 1, Ihe upper limil on Ihe cross-section can be interpreted as a limil on Ihe coupling conslanls 
divided by Ihe scale of new physics: KugtIA < 5.8 x 10~^TeV“^ assuming KcgtIA = 0, and KcgtjA < 
13 X 10“^TeV“^ assuming KugtIA - 0. Dislribulions of Ihe upper limils on Ihe coupling conslanls for 
combinations of cgt and ugt channels are shown in Fig. 10(a). 

Limils on Ihe coupling conslanls can also be interpreted as limils on Ihe branching fractions using S{t 
qg) = c{KqgtlJ^ , where C is calculated al NLO [91]. Upper limils on Ihe branching fractions !B{t 
ug) < 4.0 X 10“^, assuming !B{t eg) = 0 and S{t eg) < 20 x 10“^, assuming !B{t ug) = 0, are 
derived and presented in Fig. 10(b). 
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Figure 9: (a) Neural-network output distribution in the signal region and (b) in the signal region with neural network 
output above 0.1. In both figures the signal contribution scaled to the observed upper limit is shown. The hatched 
band indicates the total posterior uncertainty as obtained from the limit calculation. 
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Figure 10: (a) Upper limit on the coupling constants Kugt and Kcgt and (b) on the branching fractions !B(t —> ug) and 
S{t eg). The shaded band shows the one standard deviation variation of the expected limit. 
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Source Expected 95 % CL upper limit 

[pb] 

Change in the upper limit 
[%] 

Normalisation & MC statistics 

1.5 

- 

Multi-jets normalisation and modelling 

1.8 

25 

Luminosity 

1.5 

5 

Lepton identification 

1.5 

3 

Electron energy scale 

1.6 

8 

Electron energy resolution 

1.5 

4 

Muon momentum scale 

1.5 

1 

Muon momentum resolution 

1.5 

5 

Jet energy scale 

1.6 

8 

Jet energy resolution 

1.9 

32 

Jet reconstruction efficiency 

1.5 

4 

Jet vertex fraction scale 

1.5 

3 

^-tagging efficiency 

1.5 

3 

c-tagging efficiency 

1.5 

4 

Mistag acceptance 

1.5 

2 

modelling 

1.9 

34 

PDL 

1.5 

5 

Scale variations 

1.5 

2 

MC generator (NLO subtraction method) 

1.6 

8 

Parton shower modelling 

1.5 

5 

All systematic uncertainties 

2.9 

- 


Table 3: The effect of a single systematic uncertainty in addition to the cross-section normalisation and MC statist¬ 
ical uncertainties alone (top row) on the expected 95 % CL upper limits on the anomalous FCNC single top-quark 
production —> f —> Mv. The relative change quoted in the third column is with respect to the expected limit with 
normalisation and MC statistical uncertainties only. 
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8 Conclusion 


A search for anomalous single top-quark production via strong flavour-changing neutral currents in pp 
collisions at the LHC is performed. Data collected by the ATLAS experiment in 2012 at a centre-of-mass 
energy ^fs - 8TeV, and corresponding to an integrated luminosity of 20.3 fb“^ are used. Candidate 
events for which a u- or c-quark interacts with a gluon to produce a single top quark are selected. To 
discriminate between signal and background processes, a multivariate technique using a neural network is 
applied. The final statistical analysis is performed using a frequentist technique. As no signal is seen in the 
observed output distribution, an upper limit on the production cross-section is set. The expected 95 % CL 
limit on the production cross-section multiplied by the t bW branching fraction is ^ bW) < 

2.9 pb and the observed 95 % CL limit is o-qg->t x S(t Wb) < 3.4 pb. Upper limits on the coupling 
constants divided by the scale of new physics Kugtjl^ < 5.8 X 10-3 TeV"' and < 13 X IQ-^ TeV-^ 

and on the branching fractions S{t ug) < 4.0 x lO'^ and S{t eg) < 20 x lO-^ are derived from the 
observed limit. These are the most stringent limits published to date. 
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